Gravity-driven instability in a spherical Hele-Shaw cell 

Jose A. Miranda, Fernando Parisio, Fernando Moraes 
Laboratorio de Fisica Teorica e Computacional, Departamento de Fisica, 
Universidade Federal de Pernambuco, Recife, PE 50670-901 Brazil 

Michael Widom 

Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213 USA 

(February 1, 2008) 

Abstract 

A pair of concentric spheres separated by a small gap form a spherical Hele- 
Shaw cell. In this cell an interfacial instability arises when two immiscible 
fluids flow. We derive the equation of motion for the interface perturba- 
tion amplitudes, including both pressure and gravity drivings, using a mode 
coupling approach. Linear stability analysis shows that mode growth rates 
depend upon interface perimeter and gravitational force. Mode coupling anal- 
ysis reveals the formation of fingering structures presenting a tendency toward 
finger tip-sharpening. 
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I. INTRODUCTION 



The Saffman- Taylor instability fTJ has been the object of extensive study during the last 
four decades 0. It arises at the interface separating two viscous fluids constrained to flow in 
the narrow gap between closed spaced parallel plates, a device known as Hele-Shaw cell. The 
cell thickness is smaller than any other length scale in the problem, so the flow is effectively 
two-dimensional. The instability arises either from a pressure gradient advancing the less 
viscous fluid against the more viscous one, or by gravity acting on the density difference 
between the fluids. The action of such driving-force mechanisms leads to the celebrated 
viscous fingering patterns |H§. 

Most Saffman- Taylor investigations analyze flow between flat Hele-Shaw cells. In a pre- 
vious work |3[ we started studying the Saffman- Taylor problem on curved surfaces by consid- 
ering flow in a spherical Hele-Shaw cell (figure 1). The interfacial instability was produced 
by a nonzero flow injection rate Q, and gravitational effects were completely negleted. We 
examined the effect of cell curvature on the shape of the patterns and showed that positive 
spatial curvature inhibits finger tip-splitting. 

In the present paper we focus on the influence of gravity in a spherical Hele-Shaw cell. 
We consider both Q = and Q > at fixed cell curvature. The unperturbed domain 
shape is a polar cap of some size, presenting an initially nearly circular boundary, which 
is gravitationally unstable and evolves at constant area (Q = 0), or slowly increasing area 
(Q > 0), without change of topology. 

The study of viscous flow in a nonplanar Hele-Shaw cell is of interest for both scientific 
and practical reasons. On the scientific level, the influence of spatial curvature on hydrody- 
namic flow is a matter of fundamental interest. It also provides a simple mathematical model 
to describe more general situations involving the filling of a thin cavity between two walls 
of a given shape with fluid. On the practical level, it may have applications in a number 
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of industrial, manufacturing processes, ranging through pressure moulding of molten metals 
and polymer materials PJ, and formation of coating defects in drying paint thin films ||. 

The gravity- driven instability on a spherical Hele-Shaw cell also allows one to gain in- 
sight into the properties of the dynamically similar, but more complex, Rayleigh- Taylor 
instability for flow on substrates of arbitrary shapes. One familiar example of this 

type of instability is the formation of fingering patterns when chocolate syrup drains, due to 
the action of gravity, from the top of a scoop of ice cream. Despite the apparent simplicity 
of this example, it is in fact a rather complicated three-dimensional problem ||, much less 
amenable to analytic treatment than its two-dimensional Hele-Shaw counterpart. 

The geometrically constrained spherical Hele-Shaw cell forces the flow to become essen- 
tially two-dimensional, and the interface one-dimensional. High viscosity eliminates inertial 
terms from the equations of motion, making the problem simpler yet. In contrast, the con- 
ventional Rayleigh- Taylor problem is inertially driven and three-dimensional effects become 
important 

The outline of the work is the following: section derives the nonlinear equation of 



motion including both injection and gravitational driving. Section |T| discusses the resulting 
motion. Section 111 A considers linear stability analysis for purely gravitational driving with 
Q = 0. The growth is purely exponential, and the linear growth rate depends on the 
geodesic distance from the sphere's north pole. Larger distance causes both faster growth 
and more unstable modes, though the characteristic wavelength is nearly constant. The 
nonvanishing injection case Q > is studied in section |III B| . Injection of a high viscosity 
fluid tends to stabilize the interface. For Q > linear growth is non-exponential due to 
evolution of the linear growth rate with distance. Section [ill C studies the coupling of a 
small number of modes. It is shown that for gravity-driven flow fingers have tendency to 
finger tip-sharpening. 
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II. STABILITY ANALYSIS AND MODE COUPLING 



Consider two immiscible, incompressible, viscous fluids, flowing in a spherical Hele-Shaw 
cell of thickness b (see figure 1). The effetively two-dimensional flow takes place on the 
surface of a sphere endowed with the metric || 



where a is the radius of curvature of the sphere, < <p < 2ir denotes the azimuthal angle 
measured on the sphere and < p < ira is the geodesic distance from the sphere's north 
pole. The polar angle 9 = p/a. Denote the viscosities and densities of the upper and lower 
fluids, respectively as rji, Qi and r] 2 , g 2 - Consider the case Q\ > g 2 and rji > r/ 2 , and examine 
flow in the northern hemisphere 9 < tt/2. Between the two fluids there exists a surface 
tension a. The flows are assumed to be irrotational, except at the interface. Fluid 1 is 
injected into fluid 2 through an inlet located at the sphere's north pole, at a given flow 
rate Q, which is the area covered per unit time. Fluid 2 is simultaneously withdrawn, at 
the same rate, through an outlet placed at the south pole. The acceleration of gravity is 
constant, represented by g, and points from north to south pole. 

During the flow, the fluid-fluid interface has a perturbed shape described as p = 71 = 
R(t) + ((ip, t). The interface perturbation amplitude is represented by (((p, t), and R denotes 
the time-dependent unperturbed radius 



where Co = cos(Rq/o), and Ro is the unperturbed radius at t — 0. The unperturbed shape 
is a polar cap of geodesic radius p = R, surface area A = ina 2 sin 2 (i?/2a) and circumference 
C = 27rasin (R/a). 

We express the net perturbation ((<p,t) as a Fourier series 




(1) 




(2) 
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((<P,t)= £ ( n (t)exp(iwp) } (3) 

n=— oo 

where ( n (t) denotes the complex Fourier mode amplitudes, and n—0, ±1, ±2, ... is 
the discrete azimuthal wave number. The area of the perturbed shape is kept inde- 
pendent of the perturbation ( by expressing the zeroth Fourier mode in (^) as (o(t) = 
(-l/2a) cotOR/a) £ |C„(t)| 2 . 

We consider a generalized version of Darcy's law (l|J^], adjusted to describe flow between 
concentric spheres 



b 2 



V; 



Vpj - flf sin ( - ) p 



(4) 



where v, = Vj(p,ip) and Pj = pj(p,ip) are, respectively, the velocity and pressure in fluids 
j = 1 and 2. The gradient in equation (|J) is associated with the metric (|l|) and is obtained 
from the corresponding three-dimensional expression for the gradient in spherical coordinates 
(r, 9, <£>), by keeping r = a and noting that 9 = p/a. In contrast to gravity-driven flows in flat 
Hele-Shaw cells, the gravity term in (f|) is not constant, but depends on the radial distance 
p. This is a manifestation of the cell spatial curvature together with its embedding in three 
dimensional space. 

We can exploit the irrotational flow condition to define the velocity potential Vj = — V(j)j. 
Using the velocity potential, we evaluate equation (|4]) for each of the fluids on the interface, 
subtract the resulting equations from each other, and divide by the sum of the two fluids' 
viscosities to get 

A\ I -I I = -a (K)|^-7acos^— J , (5) 

where A=(r] 2 — ^i)/(?72 + ^i) is the viscosity contrast, a = b 2 a/[12{rji + r] 2 )} contains the 
surface tension, and 7 = b 2 g (qi — ^ 2 )/ [12 (771 + r/2)] is a measure of gravitational force. To 
obtain (^) we used the pressure boundary condition p 2 — p\ = <tk at the interface p = 71, 
where k is the interfacial curvature [Hi. 



Following steps similar to those performed in Ref. ||, we define Fourier expansions for 
the velocity potentials, which obey Laplace's equation. We express <pj in terms of the 
perturbation amplitudes ( n by considering the kinematic boundary condition for flow on 
a sphere. As in the flat cell case, this condition refers to the continuity of the normal 
velocity across the fluid-fluid interface [K|. Substituting these relations into equation (Q), 
and Fourier transforming, yields the mode coupling equation of the Saffman- Taylor problem 
in a spherical Hele-Shaw cell, taking into account both injection and gravity 



( n = A(n) (n + [F(n,n) Cn'Cn-n' + G(n,n) ( n '( n 



(6) 



where 



A(n) 



Q 
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(7) 



is the linear growth rate, and 
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- sgn(nn')} — C} 



(9) 



are the second-order mode coupling terms, with S = sm(R/a) and C = cos(R/a). The 
overdot denotes total time derivative and the sign function sgn(nn') = 1 if (nn') > and 
sgn{nn' 



^ - -1 if (nn') < 0. 



III. DISCUSSION 



A. Linear stability analysis with Q = 



First concentrate on the purely gravity-driven case, Q = 0. In this situation, we consider 
flow in a closed cell obeying mass conservation. To drive the interface gravitationally, one 
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could first allow the system to form a stable, unperturbed spherical cap at the south pole, 
and then invert the spherical cell to put the denser fluid on top in the unstable position at 
the north pole. 

We begin by investigating the dispersion relation (|7]). Notice that the gravity term 
in Eq. (0) contains no explicit dependence on radial distance. This finding is somewhat 
surprising, since the gravity term in Darcy's law (|J) clearly presents a dependence on p . 
Physical intuition suggests a factor of S should multiply the gravity term in (|7]), making 
it vanish at the poles and become maximal at the equator. This apparent missing factor 
reappears if we rewrite the linear growth rate in terms of the wave number k = \n\/aS. 
The variables n and k are both useful: for example, n occurs in integer values and counts 
the number of fingers. On the other hand, k determines the characteristic wavelength of a 
perturbation. 

To better understand the physical information behind the description of the linear stage 
in terms of n and k, we plot the linear growth rate at a sequence of radial distances R = a6 
in figures 2a and 2b. We use typical experimental parameters given in a recent experimental 



work in rotating, flat Hele-Shaw cells jy]]: fluid 1 is a silicone oil (r/i ~ 0.5 g/cm s, £>i ~ 1.0 
g/cm 3 ) and fluid 2 is air (r/ 2 « 0, Qi ~ 0). The thickness of the cell b = 0.1 cm and the 
surface tension a = 20.7 dyne/cm. We set the radius of the sphere a = 5 cm and acceleration 
of gravity g = 980 cm/s 2 . 

Comparing figures 2a and 2b, it is clear that the fastest growing mode (number of fingers 
n developing at a maximum growth rate) moves to large n for large radial distance (and 
small n for small radial distance) in order to keep the corresponding fastest growing k, and 
thus the apparent wavelength, nearly constant. The same effect occurs in radial flow in flat 



space [12| . In addition to the shift to larger n, the growth rate of the fastest growing mode 
increases with radial distance because the net gravitational force depends on R. This differs 



7 



from the case of radial flow in flat space IJ12 ], where the interface velocity falls off as 1/R 



causing the growth rate of the fastest growing mode to decrease. 

Consider the purely linear contribution, which appears as the first term on the right hand 
side of equation @. The condition Q = simplifies the theoretical description: R(t) = R = 
constant, and consequently the linear growth rate X(n) is time independent. This implies 
that the actual relaxation or growth of mode n is purely exponential 

Cf = 0) (t) = C»(0) exp[A(n) t], (10) 



where C n (0) is the initial perturbation amplitude. To see the overall effect of Eq. (|T0|), we 
plot evolved interfaces using the same experimental parameters as those used in figures 2a 
and 2b. It is convenient to rewrite the net perturbation (^) in terms of cosine and sine modes 
C(0,t) = Co + E [a>n(t) cos(n8) + b n (t) sin(n6 1 )], where a n = ( n + (_ n and b n = i{( n - C-n) 

n>0 

are real- valued. We take into account modes n ranging from n = 1 up to 20. Figure 3 
depicts the evolution of the interfaces, for a random choice of phases, at time t = 5 s. We 
evolve from two distinct initial radii (a) Rq = a 7r/8, and (b) Rq = a 7r/4. In both cases 
ICn(0)| = 0.05 cm and we use the same randomly chosen phases. It is evident that for larger 
radial distances (or equivalently, larger polar angles) we have both faster growth and more 
unstable modes, though the characteristic wavelength is not very different, in agreement 
with the predictions made from figures 2a and 2b. 



B. Linear stability analysis with Q > 

Now consider the nonvanishing injection case. Taking Q > introduces two important 
new effects. First, inspecting equation (|7|), we see that Q multiplies a term linear in n that 
must be added to the growth rate A obtained with Q = 0. If the inner fluid is high viscosity, 
so that A < 0, then this new term in A is negative. All modes grow more slowly with Q > 0, 



8 



and some may even become stable. Thus, Q > diminishes the strength of the instability. 
Essentially, the flow fills in the gaps between the fingers. 

The second effect is more subtle. Since the unperturbed interface radius R(t) is now 
time dependent, the linear growth rate A evolves with time as the radial distance steadily 
grows. At any instant the interface evolves exponentially with an instantaneous growth rate 
A depending on the current R(t). For sufficiently small Q the instantaneous A at R(t) nearly 
equals the Q = value for the same R(t). We call this a quasistatic approximation, and 
within this approximation it is clear that an increasing number of modes become unstable 
as time progresses due to the steady increase of R(t). The evolution of A in the quasistatic 
approximation is given by the series of curves shown in figure 2. 

If we wish to study the linear growth without making the quasistatic approximation, we 
may integrate the equation of motion (H) exactly, keeping only terms of first order in ( on 
the right-hand side. Because A is time dependent, through the variation of R(t), the exact 
solution of the linear equation of motion becomes non-exponential and can be written as 



CP >0) (t) = Cn(0) f 



tan(#/2a) 
tan(i? /2a) 



A\n\ 



exp 



27r|n| 



a{n 2 - 1) (C C, 



^r-^r) ' -C) 



(11) 

where Sq = sin (R Q /a), and R, C and S all depend on time. 



C. Mode coupling analysis 

We use the mode coupling equation (|6|) to investigate the effect of the nonlinear terms 
in the interface evolution. We are interested in studying how gravity influences the shape 
of the fingering structures, focusing on finger-tip behavior. Finger tip-splitting and tip- 
sharpening phenomena are related to the influence of a fundamental mode n on the growth 
of its harmonic In. Without loss of generality we may choose the phase of the fundamental 
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mode so that a n > and b n = 0. We replace the time derivative terms a n and b n by A(n) a n 
and A(n) 6 n , respectively, for consistent second order expressions. Under these circumstances 
the equation of motion for the harmonic cosine mode becomes 

a 2n = A(2n) a 2n + - T(2n, n) a\ (12) 



where the finger-tip function is 
T(2n,n) 



aS 1 27ra 2 S 2 



4n 



, SaC (2n 2 - 1) ) 

+ 2^ r 1 6) 



The corresponding equation for the sine modes b 2n is not as interesting as (0), since growth 
of &2n is uninfluenced by a n . 

The sign of T(2n,n) dictates whether finger tip-splitting or finger tip-sharpening is fa- 
vored by the dynamics [fL2]| . If T(2n,n) < 0, at second order the result is a driving term of 
order a 2 forcing growth of a 2n < 0. With this particular phase of the harmonic forced by 
the dynamics, the n outwards-pointing fingers of the fundamental mode n tend to split. In 
contrast, if T(2n, n) > growth of a 2n > would be favored, leading to outwards-pointing 
finger tip-sharpening. 

A noteworthy point about equation fljjD is that it shows no dependence whatsoever on 
gravity. In the evaluation of T(2n, n) from equation @ we found that the term involving 
gravity in F(2n,n) exactly cancels against the term involving gravity in X(n)G(2n,n). The 
second order term driving tip-splitting in equation (|L^) is therefore independent of the force 
of gravity, though gravity does generally influence mode coupling at second and higher 
orders. 

Inspecting Eq. (|T3|) we find that, since A = — 1 and C > 0, the finger-tip function 
T{2n,n) > for both Q = and Q > 0. Eq. (p~3|) predicts that gravity-driven flow 
on a sphere leads to patterns showing enhanced finger tip-narrowing. Informal studies of 
chocolate syrup fingers on ice cream scoops support this claim. 
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FIGURE CAPTIONS 



FIG. 1: Schematic configuration of flow in a spherical Hele-Shaw cell. 

FIG. 2: Linear growth rate as a function of (a) mode number n and (b) wave number k 
at a sequence of radial distances (1) R — a 7r/8, (2) R — a ir/4, (3) R — a 3ii/8 and (4) 
R = a tt/2. 

FIG. 3: Interface evolution according to Eq. (|1(]) at t — 5 s, for (a) R = a tt/8 and (b) 
Ro = a 7r/4. Other parameters are given in the text. 
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